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Abstract 

Neutrino masses and mixings are analyzed in terms of left-handed 
fields and a 6 x 6 complex symmetric mass matrix Ai whose 
singular values are the neutrino masses. An angle character- 
izes the kind of the neutrinos, with = for Dirac neutrinos 
and = tt/2 for Majorana neutrinos. If Xjy = 0, then baryon- 
minus-lepton number is conserved. When x^, ~ 0, the six neu- 
trino masses coalesce into three nearly degenerate pairs. Thus 
the smallness of the differences in neutrino masses exhibited in 
the solar and atmospheric neutrino experiments and the stringent 
limits on neutrinoless double beta decay are naturally explained if 
B — L is approximately conserved and neutrinos are nearly Dirac 
fermions. If one sets sinx^, = 0.003, suppresses inter-generational 
mixing, and imposes a quark-like mass hierarchy, then one may 
fit the essential features of the solar, reactor, and atmospheric 
neutrino data with otherwise random mass matrices Ai in the 
eV range. 

This B — L model leads to these predictions: neutrinos oscillate 
mainly between flavor eigenfields and sterile eigenfields, and so 
the probabilities of the appearance of neutrinos or antineutrinos 
are very small; neutrinos may well be of cosmological importance; 
in principle the disappearance of Vr should be observable; and 
Oi//3/3 decay is suppressed by an extra factor of 10~^ and hence 
will not be seen in the Heidelberg/Moscow, IGEX, GENIUS, or 
CUORE experiments. 
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1 Introduction 



There are three major sets of experimental results that shed light on neutrino 
masses and mixings. Solar-neutrino and reactor experiments have shown that 
the electron neutrino z/g couples to two nearly degenerate mass eigenfields 
with 10~^° ^ ml — m\ <, 10^^ eV^. The atmospheric-neutrino experiments 
have shown that the muon neutrino and antineutrino P^^ couple to two 
nearly degenerate mass eigenfields with 10~^ ^rnl — ml^ 10~^ eV^. Double- 
beta-decay experiments have placed very stringent limits on the rates of 
neutrinoless double beta decay. 

The current wisdom on neutrinos is that they have very small masses 
and that the tiny mass differences experimentally observed are due to their 
small masses. It is also generally believed that the smallness of the masses 
of the neutrinos arises from the seesaw mechanism [1], which involves huge 
Majorana masses. It is also thought that CKM-hke mixings explain the solar 
and atmospheric data. 

The purpose of this paper is to present a general discussion of neutrino 
masses and mixings and a rather different explanation of the experimental 
facts. It is argued that the neutrino mass matrix is a 6 x 6 complex symmetric 
matrix, that the moduli of the Majorana mass terms are small compared to 
those of the Dirac mass terms, that neutrinos propagate as six mass eigen- 
fields with masses that form three nearly degenerate pairs, and that the 
oscillations observed in the solar and atmospheric experiments are mainly 
between members of these pairs. Neutrino masses and mixings have been 
discussed by several authors [2, 3], and the model introduced by Geiser [4] is 
an important example of the class of models to be developed in what follows. 

Because left-handed fields can mix only with other left-handed fields, 
neutrino masses are best analyzed in terms of left-handed fields. There are 
six such fields, the three left-handed fiavor eigenfields t/g, and Vr-, and the 
charge conjugates of the three right-handed neutrinos, which are sterile. In 
section 2 the neutrino mass matrix is exhibited as a 6 x 6 complex symmetric 
matrix M.. In general this matrix is not normal, but like every matrix it 
admits a singular- value decomposition M. = U MV'^ in which the matrices U 
and V arc unitary and the matrix M is real, nonnegative, and diagonal [5]. 
Moreover because the matrix }A is symmetric it follows from a theorem of 
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Takagi's [6, 7] that one may always choose V = U*. In section 3, Takagi's 
factorization Ai = UMU^ is used to diagonaUze the action and to show 
that the singular values of the neutrino mass matrix Ai or equivalently the 
diagonal elements of the matrix M are the six neutrino masses mj and that 
the unitary matrix U describes the neutrino mixings. 

In section 4 the LEP measurement of the number of light neutrino fla- 
vors is related to the mixing matrix U. Section 5 describes a cosmological 
constraint on the masses of light, stable neutrinos. Section 6 expresses the 
neutrino oscillation probabilities in terms of the mixing matrix U and the 
masses rrij and briefly discusses the results of the solar, reactor, and atmo- 
spheric neutrino experiments. 

In section 7 an angle Xj, is introduced that quantifles the extent to which 
neutrinos are Dirac fermions or Majorana fermions. Dirac neutrinos have 
Xi, — 0, and Majorana neutrinos have Xi, = ti/2. If all Majorana mass terms 
vanish, that is if x^, = 0, then the standard model conserves baryon-minus- 
lepton number, B—L, which is a global U (1) symmetry. It is therefore natural 
in the sense of 't Hooft [8] to assume that x^, ~ so that this symmetry 
is only slightly broken. The neutrinos then are nearly Dirac fermions and 
their masses coalesce into three pairs of almost degenerate masses. Thus the 
approximate conservation oi B — L explains the tiny mass differences seen 
in the solar and atmospheric neutrino experiments without requiring the 
neutrino masses to be absurdly small. There are three tiny mass differences, 
of which one explains the solar experiments and another the atmospheric 
ones, and three unconstrained mass differences, which can lie in the cV range. 
If one sets sina;,^ = 0.003, suppresses inter-generational mixing, and imposes 
a quark-like mass hierarchy, then one may fit the essential features of the 
solar, reactor, and atmospheric neutrino data with otherwise random mass 
matrices Ai in the eV range. Thus neutrinos easily can have masses that 
saturate the cosmological bound of about 8 cV. Moreover because neutrinos 
are almost Dirac fermions, neutrinoless double beta decay is suppressed by 
an extra factor ~ sin^ x^, sin^ yi, < 10~^, where yu is a second neutrino angle, 
and is very slow, with hfetimes in excess of 2 x 10^^ years. 

An appendix outlines an efficient way of performing the singular-value 
decomposition of the mass matrix M. by means of a call to the LAPACK [9] 
subroutine ZGESVD [10]. 



3 



This B — L model of neutrino masses and mixings leads to these predic- 
tions about future experiments: 

• The three flavor neutrinos oscillate mainly into the conjugates of the 
right-handed fields, which are sterile, but the mass difi:erenccs among 
the three pairs of nearly degenerate masses can lie in the eV range. The 
probabilities of the appearance of neutrinos or antineutrinos are small, 
as shown by LSND and KARMEN2, but MiniBooNE has a reasonable 
chance of seeing the appearance of ^ i/g- 

• Because neutrino masses are not required to be nearly as small as the 
solar and atmospheric mass differences might suggest, they may well 
be of cosmological significance. 

• If a suitable experiment can be designed, it should be possible to see 
the tau neutrino disappear. 

• The rate of neutrinoless double beta decay is suppressed by an extra 
factor of ~ sin^ Xi, sin^ yi, ^ 10~^ and hence will not be seen in the 
Heidelberg/Moscow, IGEX, GENIUS, or CUORE experiments. 



2 Mass Terms and Mass Matrices 

There arc three potential sources of neutrino masses. One source is Dirac 
mass terms, which require the existence of right-handed neutrino fields. The 
second source of neutrino masses is Major ana mass terms composed of left- 
handed neutrino fields; these mass terms require a triplet of Higgs bosons, 
break B — L, and drive neutrinoless double beta decay. The third source is 
Majorana mass terms composed of right-handed neutrino fields; these terms 
also break B — L but do not affect neutrinoless double beta decay, at least 
in leading order. Inasmuch as quarks and charged leptons come in both left- 
handed and right-handed fields, the existence of right-handed neutrino fields 
is reasonable. And because the Higgs sector is still unknown, the existence 
of Higgs bosons that transform as a triplet under SU (2) l remains possible. 
So in this paper we shall consider all three kinds of mass terms. 
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Because left- and right-handed fields transform differently under Lorentz 

boosts, they cannot mix. It is therefore convenient to write the action den- 
sity exclusively in terms of two-component, left-handed fields. The two- 
component, left-handed neutrino flavor eigenflelds z/g, '^/i, i^r will be denoted 
z/j, for i — e, II, r. The two-component, left-handed fields that are the charge 
conjugates of the putative right-handed neutrino fields rire, nj.^, rirr will be 
denoted rij = — io"2 n*.- for i = e, fi,T in which hermitian conjugation with- 
out transposition is denoted by an asterisk, and a2 is the second Pauli spin 
matrix. 

For two left-handed fields ip and x, there is only one mass term, and its 
hermitian form is 

2 2 

i mx'^CTat/' - i m^^/'VsX* = ^ ("^ Xa(^2api^p - m* C^2a/3X^) , (1) 

a=l (3=1 

which incidentally is invariant under a Lorentz transformation g G SL(2, C) 
because the antisymmetric matrix (72 converts two factors of g into det(5') = 
1. If the field x is the charge conjugate of a right-handed field that carries 
the same conserved quantum numbers as the field ip, then this mass term (1) 
is called a Dirac mass term. On the other hand, if the fields ip and x carry 
the same quantum numbers or no conserved quantum numbers at all, then 
the mass term (1) is called a Majorana mass term. Since electric charge is 
conserved, only neutral fields can have Majorana mass terms. 

The six left-handed neutrino fields UijUi for i = 1,2,3 can have three 
kinds of mass terms. The fields Ui and rij can form Dirac mass terms 

iDijuJa2nj - iD*p]a2V* (2) 

in which the Dij are complex numbers. In a minimal extension of the stan- 
dard model, the Dij are proportional to the mean value in the vacuum of 
the neutral component of the Higgs field. The fields rii and rij can form 
Majorana mass terms 

iEijnJa2nj - iE*jn\a2n*. (3) 

Within the standard model, the complex numbers E^j are simply numbers; 
in a more unified theory, they might be the mean values in the vacuum of 
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neutral components of Higgs bosons. The fields i/j and uj can also form 
Majorana mass terms 

iF,ji^Ja2i^j-iF*jU^a2U*. (4) 

In a minimal extension of the standard model, the complex numbers Fij might 
be proportional to the mean values in the vacuum of the neutral component 
of a new Higgs triplet hab = hha that transforms as 

Khi.^) = K'b'ga'ii.x)gvl{x) (5) 

in which the 2x2 matrix Qaa' {x) is the one that transforms the three doublets 

^^(^^^(eiS) ^ Hx)' ^ g{x)U{x) (6) 

where Cj = e, ii,t for i — 1,2,3. Such a triplet Higgs multiplet would allow 
a term like 

2 

Yl hab{x)Lia{xy a2Ljb{x) (7) 

a,b=l 

to remain invariant under SUii^) ® U{1)y gauge transformations and so to 
contribute a mass term to the action density. 

Since (J2 is antisymmetric and since any two fermion fields x and ip anti- 
commute, it follows that 

x'^(72'0 = ip^cr2X and xV2'0* = ip''cr2X*: (8) 

which incidentally shows that a singlet Higgs field, hab = —hba, would lead to 
Fij — 0. This symmetry relation (8) implies that the 3x3 complex matrices 
E and F are symmetric 

E'^^E and F^^F (9) 

and that 

iDijnJa2iyi = iDijvJ a2nj. (10) 
Let us introduce the 6x6 complex matrix 

M = [^, b) (11) 
6 



which by (9) is symmetric. Then with the six- vector N of left-handed neu- 
trino fields 



N 



He 

n„ 



(12) 



we may gather the mass terms into the matrix expression 



-N^Ma2N 
2 



-N''M*a2N*. 



(13) 



Actually since each of the six fields in TV has two components, the six-vector 
N has twelve components. The use of twelve-component fields has been 
advocated by Rosen [11] and more recently by Starkman and Stojkovic [12]. 

The complex symmetric mass matrix }A is not normal unless the positive 
hermitian matrix A4A4^ is real because the commutator [M., M'^] is twice its 
imaginary part: 

[M, A^^] = 2i (MM^) . (14) 

When the mass matrix is itself real, then it is also normal and may be 
diagonalized by an orthogonal transformation. But in general A4 is neither 
real nor normal. 

Yet like every matrix, the mass matrix M. possesses a singular- value de- 
composition [5, 9] 

M = UMV^ (15) 

in which the 6 x 6 matrices U and V are unitary and the 6x6 matrix M is 
real and diagonal 
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(16) 



with singular values rrij > 0. These singular values will turn out to be the 
masses of the six neutrinos. The columns of the unitary matrix U are the 
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left singular vectors of M.; they are the eigenvectors of the matrix M.A4^ 
with eigenvalues m| 

MM^U M\ (17) 

The columns of the unitary matrix V arc the right singular vectors of Ai; 
they are the eigenvectors of the matrix M.^ A4 with eigenvalues m| 

M'' MV M^. (18) 



We have not yet made use of the fact that the complex matrix M. is 
symmetric. It is a general mathematical theorem [6] due to Takagi [7] that 
for every symmetric complex matrix Z there is a unitary matrix X and a real 
diagonal matrix D with nonnegative elements such that Z = XDX^ . The 
diagonal matrix elements Djj are the singular values of the matrix Z. Thus 
the symmetric complex matrix A4 may be factored by a unitary matrix W 

M = WMW^ (19) 

where M is the matrix of masses (16). The columns of W are the eigenvectors 
of the matrix A4 A4^ with eigenvalues ntj 

MM'^W = W M^. (20) 

Thus by (17) the columns of W are the same as the columns of the unitary 
matrix U of the singular-value decomposition (15) apart from over-all phase 
factors. Similarly the columns of W* are the eigenvectors of the matrix 
Ad'^ M. with eigenvalues m| 

MW* ^W* M'^, (21) 

and so by (18) the columns of W* are the same as the columns of the uni- 
tary matrix V of the singular-value decomposition (15) apart from over-all 
phases. For our purposes the import of Takagi's theorem is that the singular- 
value decomposition of the symmetric complex matrix Ai — U MV^ may be 
achieved with U and V = W* so that 

M^UMU^, (22) 

which we may call Takagi's factorization. 
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3 Field Equations, Masses, and Mixings 



The free, kinetic action density of a two-component left-handed spinor t/j is 
iip'^ {do — a ■ V) ip. Thus by including the mass terms (13), one may write 
the free action density of the six left-handed neutrino fields N as 

Co = iN^ {do-a-V)N+ -N^Ma2N - -N^M*a2N* (23) 
or with Takagi's factorization (22) of the mass matrix as 

Co = iN^ {do-a-V)N+ ^N^UMa2 N - ^N^U* Ma2U^ N* (24) 

Zi Zi 

in which M is the real 6x6 diagonal matrix (16) with nonnegative entries 
rrij. Let us define the six- vector of fields 



Nm^U^N or Vmi^Y^UjiNj. (25) 

Since the 6x6 matrix U is unitary, we have 

N = U*Nm and = N^U^ (26) 

as well as 

N*^UN:^ and N^^NlU^. (27) 
Thus we may write Co in the form 

Co = iNl {do- a- V) U* 

+ l-Nl UM(J2U^ U* - l-Nl U*M(T2U^ U AT;; (28) 

Zi z 



or since U'^ U* ^ I ^ UW 

Co = iNl {do- a- V) + '-Nl Ma2 - '-NI Ma2 N^. 
The action density therefore has the diagonal form 

6 r ^ 

^0 = ^ i^l^ {do- a- V) u^^ + -rrii (i/^. as i^m^ - i^L, (^2 i^^,) 
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i=l 



(29) 



(30) 



The fields Ujn- are the normal modes of the theory. They propagate ac- 
cording to the equation of motion 

{do- a- V) = rUi as (31) 

in which the asterisk represents complex or hermitian conjugation without 
transposition. By taking the hermitian adjoint of this equation without trans- 
posing the matrices and vectors (or equivalently by taking the hermitian ad- 
joint of the equation and then re-transposing the matrices and vectors) , one 
has 

{do-a*-W)u*^^^m,a;u^,. (32) 

The rules a2a*a2 = — c?, (c2)^ = 1, and a^* — —ui now imply that the 
equation of motion for i/j^. is 

(9o (J • V) (yiv*m^ = -Vdi Urrn ■ (33) 

Applying {do + a ■ V) to the field equation (31) for i/^. and then using the 
field equation (33) for u^., we find that 

(□ - m^) vm, = 0. (34) 

Thus the diagonal elements rrii of the real mass matrix M, which are the 
singular values of the mass matrix are the neutrino masses, and the 
eigenfield of mass is as given by (25). A neutrino of mass described 
by the two-component field must have two helicities or spin states and so 
must be its own antiparticle. The two-component field is not hermitian, 
however; we'd need to use four components to make this field hermitian. 

Since the mass eigenfields Vrrn are the normal modes of the theory, it is 
worth recoding their expansions 



^mi{x) = j J^^^J^ ^(P' nii)e'^''a{p, s, rrii) + v{p, s, mi)e '^''a\p, s, 

(35) 



(27r)3/2 

in which px = p- x — et. In terms of the normalization factor 



m = (36) 
2^e{e + mi) 
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the two-component spinors u and v are 

^p,i,m,)=n,(^_: + ^T^^^) u{p,-lm,)=n,^-Pl+^'J (37) 
and 

v{p,^,mi) ^ rii ( \ ^(-^ 1 j ^ / P3 e\ ^^^^ 

They satisfy the spin sums 

u{p, s, rrii) u\p, s, rui) = ^ v{p, s, rrii) v^p, s, m^) = ^{e - p- a) (39) 

s s 

and 

Yu{p, s, mi)v^{p, s, rrii) = (J2. (40) 
s 2e 

For momentum p = pz in the z- direction with p ^ rrii, these spinors 
reduce to 

»(p.i."..)-|(i) 4. «.)«(;) (41) 

and 

v{p,\,mi) ^ ^;(p,-i,mi) ~ , (42) 

which shows that the field I'mii^) primarity deletes neutrinos of negative 
helicity and adds neutrinos of positive helicity. 

In terms of the flavor eigenflelds i/j and n^, the eigenfield of mass rrii is 

6 3 3 

^ruj = X] ^a^i = Yj ^ij^i + U^i+3)jni. (43) 

i=l i=l 1=1 

The flavor eigenflelds N are given by 

N - U*Nm. (44) 

In particular, the three left-handed fields Ui for i = e, ij,,t are linear combi- 
nations of the six mass eigenfields 



-Y.U:,yra, (45) 
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and not simply linear combinations of three mass eigenfields. The three 
active (right-handed) antineutrino fields u* for i — e, ii,t are similarly 

6 

ia2iy* = Y,Uijia2iy;^.. (46) 

The matrix that expresses the flavor eigenflelds in terms of the mass eigen- 
fields is a 3 X 6 matrix that is half of a 6 x 6 unitary matrix, not a 3 x 3 
unitary matrix. The three sterile left-handed fields rii for i — e, ii,t are also 
linear combinations of the six mass eigenfields u^. 

6 

^i = 12^+3)j^rnj- (47) 



In view of the relation (45) between the flavor eigenflelds and the mass 
eigenflelds, the action density for neutrino interactions in four-component 
form [13] 



4 sin ~r[ 
may be written as 



^ a/2 [e, r (1 + l^>i + + 75)e,] - ^ 



Vi ^(1 + 75)z/i 



COS 9 



w 



(48) 



-le 



4 sin 9w ~y 



6 6 

ei ^(1 + 75) E ^''^i + E ^ik^mu + 75)ei 

i=\ k=\ 



j=l fe=l 



COS 9 



w 



(49) 



In two-component form, this is 



2sm9w^,^^ 



i=i 



k=l 
6 6 



cos 9\ 



j=lk=\ ^^^JW 

in which all flelds are left-handed two-component spinors. 



(50) 
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4 LEP 



The four LEP measurements of the invisible partial width of the Z impose 
upon the number of light neutrino types the constraint [14] 

= 2.984 ± 0.008. (51) 

It follows from the action density (49) that the amphtude for the Z pro- 
duction of two neutrinos v^^^ and i/^^. to lowest order is 

3 

A{v^^,Vra^) (xY.UikU*j, (52) 
i=l 

and therefore that the cross-section for that process is 

3 

0-(i^m,-,i^mJ OC I ^C/ifeC/* |^. (53) 
1=1 

The measurement (51) of the number Ni, of light neutrino species thus implies 
that the sum over the light-mass eigenfields is 

3 

E I E^^fe^r/ = 2.984 ±0.008. (54) 

j,k light 1=1 

This constraint on the 6x6 unitary matrix V is quite well satisfied if all 
six neutrino masses are light. For in this all-light scenario, the sum is 

E E E U,kU;pkUi'j = jljl ^w^vv = E l = 3 2.984 ± 0.008. (55) 

j,k=l i=l i'=l 1=1 i'=l i=l 

But the constraint (51) will also be satisfied if the three flavor eigenfields 
couple only to light-mass eigenfields z/^ , for in this case the matrix elements 
U*j between flavor i and heavy mass mj vanish, and one may extend the sum 
over the light-mass eigenfields to a sum over all six mass eigenfields. In this 
few-light scenario, the three independent fiavor eigenfields Ui must couple to 
at least three hght-mass eigenfields . . 
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5 Cosmology 



5.1 Weighing Neutrinos 

Any stable, two-component neutrino mass eigenfield that couples via 
U*j to a flavor eigenfield must have a contemporary density of 115/cc due 
to thermal equilibrium in the early universe [15, 16]. Thus the present 
mass density of the eigenfield is 115mj/cc. The current critical den- 
sity is 1.05 X 10^/i^eV/cc, where h is the Hubble parameter in units of 100 
km/sec/Mpc. Thus ii h 0.65, then the contribution of to the neutrino 
critical density is rrij/ 40 eV. So the contribution to of the stable, light 
mass eigenfields that couple to the three fiavor neutrinos is 

I Y mA (56) 

Since O,^ is the neutrino contribution to Hot Dark Matter, and since the 
HDM part of Qm is probably less than 0.2, we may conservatively conclude 
that Qj, ^ 0.2. Thus we arrive at an approximate upper bound [15, 16] for 
the sum of the masses rrij of the light, stable neutrinos that couple to the 
three flavor eigenflelds: 

J2 mj < 8eV. (57) 

j light 

Under other assumptions Fukugita, Liu, and Sugiyama [17] have derived 
limits on the sum of the light neutrino masses that are in the range of 2 to 5 
eV. The Sloan Digital Sky Survey [18] will measure fli, and weigh the light 
neutrinos [19] . Within the context of a minimally extended standard model, 
these cosmological bounds and the LEP constraint (54) imply that the Z 
gauge boson couples to between 3 and 6 very light neutrinos. 



5.2 BBN Bounds on Light Sterile Neutrinos 

It is open question whether active neutrinos oscillate mainly into other active 
neutrinos or mainly into sterile neutrinos. If they oscillate mainly into sterile 
neutrinos, then in the early universe electron neutrinos before they decoupled 
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might have brought sterile neutrinos into chemical equilibrium [20], increas- 
ing the effective number N'^f of two-component neutrino species. A higher 
Nl^ would have raised the energy density p of the universe and therefore its 
rate R of expansion since E? = {8nG/3)pR'^ — k where k is the curvature. 
The neutron-to-proton ratio would then have been higher at its freeze-out 
(T Ki 0.7 MeV), and the abundance Y of helium higher than is observed 
{Y ^ 0.24 ± 0.01). In the case of maximal mixing between z/^ and u^g, the 
conservative BBN limit A^"^^ < 4 would impose on Sml^ = \m1 — ml ^\ 
an upper limit of the order of 10~^ eV^ [20]. This big-bang-nucleosynthesis 
(BBN) constraint would rule out the possibihty that the atmospheric muon 
neutrinos oscillate mainly into sterile muon neutrinos. 

5.3 Lepton Asymmetries 

However an excess of Ue over Pe, or vice versa, when the temperature was 
between 30 and 0.7 MeV, would have modified the reaction rates n -|- i^e <-> 
p + e~ and n + e'^ ^ p + and consequently changed the effective number 
of neutrino species N^^ . Thus Foot and Volkas have pointed out [21] that 
the BBN constraints on active-sterile neutrino oscillations depend upon the 
assumption tliat the lepton-number asymmetry of the early universe was 
neghgible. They also have argued [21] that active-sterile neutrino oscillations 
could themselves have generated a significant lepton-number asymmetry at 
the relevant temperatures. So the BBN limits on N^^ and Srn^^ may not 
be valid. And the recent measurement of the cosmic microwave background 
(CMB) radiation by the Boomerang experiment has called into question the 
precision of BBN results; it now appears that the BBN bound on the baryon 
density ought to be raised by about 50% to 0.024 < h^Q^ < 0.042 [22]. 

The Affleck-Dine mechanism [25] has been suggested as a second way of 
generating B and L asymmetries in the minimal supersymmetric model [27] 
and its extensions [26]. 
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5.4 Baryon Number of the Universe 

K. Dick, M. Lindner, M. Ratz, and D. Wright [23] recently showed how a uni- 
verse that conserves B — L might have started with B = L = and still have 
arrived at the observed current baryon-to-photon ratio as long as the masses 
of the neutrinos are suitably small. In their model a density oi B + L pro- 
duced at GUT temperatures is equilibrated between left- and right-handed 
(then massless) particles by Yukawa processes at temperatures above 1 TeV. 
Sphalerons [24] then washed out left-handed baryons and leptons, reducing 
B + L while conserving B — L. If the Yukawa couplings of the neutrinos 
are so weak as to yield (purely Dirac) neutrino masses less than about 10 
keV, then the equilibration of left- and right-handed neutrinos could not have 
kept up with the sphaleron washout, and an excess of left-handed baryons 
equal to the excess of right-handed neutrinos, which are immune to sphaleron 
washout, would have survived the electroweak era. This mechanism is both a 
potential explanation of the baryon number of the universe and a third way 
of generating lepton asymmetries. 

Since B — L is exactly conserved in this model, neutrino oscillations occur 
only among the three active flavors. But this model might work even \i B — L 
were slightly broken. 



6 Oscillations 

It follows from the action density (49) and from arguments presented else- 
where [29] that the lowest-order amplitude A{vi — > Vii) for a neutrino Vi {e.g., 
produced by a charged lepton Cj) after propagating with energy E a distance 
L as some light-mass eigenfleld of mass ruj <^ E to appear as a neutrino u'- 
{e.g., producing a charged lepton e[) is 

A{u^ ^ ly,) = E t/;,.t/,,e--?^/(2^). (58) 

j light 

In view of (46), the lowest-order amplitude for the anti- process, Ui — > Ui', is 
A{Pi ^ P,0 = E C/i',C/* e-^"^?^/(^^), (59) 

j light 
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which is A(z/j — > i/j/) with the matrix elements of U replaced by their complex 
conjugates [28]. 



The corresponding probabilities to lowest order are 



and 



P{i^i ^ = Yl U;jUi,UefU*j, exp [i{m} - m])L/{2E) 

light 



P(j^. ^ i,.,) = ^ U;jUijUi,fU*j, exp [z(m2 - m^j,)L/{2E) 

light 



(60) 



(61) 



both of which, if all six neutrinos are light, approach 5^ in the limit L/£' — > 0. 

A measure oi CP violation is provided by the asymmetry parameter 
A{i,i') defined as 



P{Ui i/j/) - P{Ui i/j/) 



and given by 



light 


Im 




1 sin 


"(m| - m})L/{2E) 


j' light 


Re 


[u*,,u,,u,,.u;,] 


cos 


\m] - m})L/{2E) 



(62) 



(63) 



As we shall see in the next section, if the six neutrinos are light and either 
purely Dirac or purely Majorana, then for each i = e, ij,,t the sum over i' is 
unity: 

j2P{u,^u,,) = l; (64) 

i'=e 

but if the neutrinos arc nearly but not quite Dirac fermions, then this sum 
of probabilities tends to be about 



1 



(65) 



If for simplicity we stretch the error bars on the Chlorine experiment, 
then the solar neutrino experiments, especially Gallex and SAGE, show a 
diminution of electron neutrinos by a factor of about one-half: 



(66) 
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which requires a pair of mass eigenstates whose squared masses differ by at 
least ~ 10^^° cV^ [15]. The reactor experiments, Palo Verde and especially 
Chooz, imply that these squared masses differ by less than ~ 10~^eV^ [15]. 

The atmospheric neutrino experiments, Soudan II, Kamiokande III, IMB- 
3, and especially SuperKamiokande, show a diminution of muon neutrinos 
and antineutrinos by about one-third: 

2 

Patm(«^/. ^ Z^m) ~ g, (67) 

which requires a pair of mass eigenstates whose squared masses differ by 
10-3 < \m] -ml\< 10-2 eV^ [15]. 

If the LSND neutrino experiment is correct, then it requires a pair of 
states whose squared masses differ by 10~^ ^Irrij — ml\ < 10+^ eV^ [15]. 



7 The B- L Model 

When the Majorana mass matrices E and F are both zero, the action density 
(23) is invariant under the U{1) transformation 

jy'^^iOGjy (gg) 

in which the 6x6 matrix G is the block-diagonal matrix 

«=(o -°.) (-) 

with / the 3x3 identity matrix. The kinetic part of (23) is clearly invariant 
under this transformation. The mass terms are invariant only when the anti- 
commutator 

F 



vanishes. 



{M,G} = 2[^^ _^]=0 (70) 



This U{1) symmetry is the restriction to the neutrino sector of the sym- 
metry generated by baryon-minus-lepton number, B — L, which is exactly 
conserved in the standard model. A minimally extended standard model 
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with right-handed neutrino fields n^j and a Dirac mass matrix D but with 
no Majorana mass matrices, E = F = 0, also conserves B — L. When B — L 
is exactly conserved, i.e., when D 7^ but E = F = then the six neu- 
trino masses rrij collapse into three pairs of degenerate masses which may be 
combined into three Dirac neutrinos. 

Suppose this symmetry is slightly broken by the Majorana mass matrices 
E and F. Then for random mass matrices D, E, and F, the six neutrino 
masses rrij will form three pairs of nearly degenerate masses as long as the 
ratio 

Tr{2DW + E^E + F^F) ^ ^ 

is small. For a generic mass matrix A4, the parameter sm^ x,y lies between 
the extremes 

< sin^ < 1 (72) 

and characterizes the kind of the neutrinos. The parameter sin^ x^, is zero 
for purely Dirac neutrinos and unity for purely Majorana neutrinos. 

Let us now recall 't Hooft's definition [8] of naturalness: It is natural to 
assume that a parameter is small if the theory becomes more symmetrical 
when the parameter vanishes. In this sense it is natural to assume that the 
parameter sin^ x^, is small because the minimally extended standard model 
becomes more symmetrical, conserving B — L, when sin^ Xi, — 0. 

In Fig. 1 the six neutrino masses rrij are plotted for a set of mass matrices 
M. that differ only in the parameter sinxi,. Apart from sinxj,, every other 
parameter of the mass matrices Ai is a complex number z = x + iy in which 
X and y were chosen randomly and uniformly on the interval [— leV, leV]. 
It is clear in the figure that when sin^Xj, ~ 0, the six neutrino masses rrij 
coalesce into three nearly degenerate pairs. Although the six masses of the 
neutrinos are in the eV range, they form three pairs with very tiny mass 
differences when sin^ ~ 0. 

Thus the very small mass differences required by the solar and atmo- 
spheric experiments are naturally explained by the assumption that the sym- 
metry generated by 5 — L is broken only slightly by the Majorana mass 
matrices E and F. This same assumption implies that neutrinos are very 
nearly Dirac fermions and hence explains the very stringent upper limits on 
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Random Neutrino Mass Matrices 




0.2 0.4 0.6 0.8 1 



sinxy 

Figure 1: The six neutrino masses are plotted against the parameter sinxj, 
for a set of random 6x6 mass matrices. 



neutrinoless double beta decay [15]. Because the masses of the six neutrinos 
may lie in the range of a few eV, instead of being squashed down to the meV 
range by the seesaw mechanism, they may contribute to hot dark matter in 
a way that is cosmologically significant. This B — L model with sin^ 
is the converse of the seesaw mechanism. 



An Example 



A simple illustration of these ideas is provided by the mass matrix A4 



M 



/ /i sin a; di cos x 

/asinx 

fssi-ax 

di cos a; 64 sin x 

ci2Cosa; 

\ ds cos X 



\ 

(^2 COS X 

ds cos X 



65 sin X 

66 sin x / 



(73) 



in which all the elements are taken to be real and non-negative. At x — 
7r/2, the Majorana parameter sinxu is unity, and the neutrinos are purely 
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Majorana. The six flavor eigenfields 1^1,1^2, 1^3, n-i,n-2,n-3 are then also mass 
eigenfields with masses /i, /2, /s, Ci, 62, 63, respectively. 



For X between 7r/2 and zero, the fields and rii for each i — 1, 2, 3 mix 
to form three pairs of mass eigenfields with masses 



+ fi) sin X ± y {ei — fi)^ sin x + Adq cos^ x 



(74) 



which are the singular values of M.. The difference of the squared masses of 
the ith pair is proportional to sin x 



— = (cj + fi) y (cj — /j)^ sin^ x + Adf cos^ x sin x. (75) 

If the mass parameters, fi,ei,di, are comparable in magnitude, then the 
Majorana parameter sin x^, is roughly sin x. 

At X — 0, Xi, — 0, the neutrinos are purely Dirac fermions, and the pair 
of mass eigenfields 

1 —i 

^mi+ ^ + ^^"^ ^-^i- ^ ^ ~ ^^^^ 

have the same mass or singular value, rui = \mi±\ = di. Since these fields are 
degenerate, the linear combinations 

+ ^^rui-) and (i/^,+ - ivrrn-) (77) 

are also mass eigenfields; they may be combined to form a single Dirac neu- 
trino 

of mass rrii for each i = 1,2,3. It follows from Eqs. (31) and (33) that ipi 
satisfies Dirac's equation {■y"'dn + mi)ipi = 0. When Xy is not zero or near 
zero, then there is no reason why any two of the masses mj should be equal, 
and one may not combine two of the fields to form a Dirac neutrino. 

For very small values of the angle x and of the parameter siv? Xy, the 
neutrinos are nearly Dirac, and the three squared-mass differences 

— rnl_ fti 2di{ei + fi) cos x sin x (79) 
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being proportional to sinx ~ sinxj,, are very small, and the mass eigenfields 
are approximately 

t^mi± ^ ^{iyi± Hi) + tanxn,. (80) 

If sin^ Xu = 0, then there are three purely Dirac neutrinos, and the mixing 
matrix V is block diagonal 

in which the 3x3 unitary matrices u and v occur in the singular- value decom- 
position of the 3x3 matrix D = umv^. If these three Dirac neutrinos are 
also light, then unitarity implies that the sum of the normalized probabilities 
is unity 

^ P{ui ^ Ui,) = 1. (82) 

i'=e 

If sin^a;,^ = 1, then this sum is also unity by unitarity because in this case 
the mixing matrix for the six purely Majorana neutrinos is 

But if there are six light, nearly Dirac neutrinos, then each neutrino flavor 
Ui will oscillate both into other neutrino flavor eigenflelds and into sterile 
neutrino eigenfields. In this case this sum tends to be about a half 

j^Piui^u,)^^. (84) 



Fitting the Solar and Atmospheric Data 

This intuition is supported by Figs. 2 and 3 which respectively display for 
10000 random mass matrices M. the sums of probabilities J2i'=e Pso\{i^e ^v) 
and Z^P=e -fatm(i^/i — ^i') ^s a function of the parameter sinx,^. The plots 
clearly show that these sums tend to cluster around the value | when sin x^, 
is small but not infinitesimal. The points at sinxj, = and at sinx,^ = 1 
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Summed Solax Probabilities for Random Mass Matrices 
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Figure 2: The sum of the probabihties Pso\{i^e ^i) for solar neutrinos 
summed over i for i = e, fi,T for 10000 random mass matrices A4 are plotted 
as a function of the parameter sinx,^. 



display the unitarity relation (82) for respectively purely Dirac and purely 
Major ana neutrinos. 

In these scatter plots and in those that follow, every parameter of the 
10000 matrices is a complex number z = x + iy with x and y chosen ran- 
domly and uniformly from the interval [— leV, leV]. The solar neutrinos are 
taken to have an energy of 1 MeV, and the probability (60) is averaged over 
one revolution of the Earth about the Sun. The atmospheric neutrinos are 
averaged over the atmosphere and over energies in the range of 1-30 GcV 
weighted by the flux of atmospheric muon neutrinos as a function of energy 
and local zenith angle given by the Bartol group in Table V of ref. [30]. 

If we continue to interpret the Chlorine experiment very conservatively, 
then the data of the solar and atmospheric experiments (66) and (67) require 
that 

1 2 

But because of the approximate, empirical sum rule (84) for ? = e and fi, the 
only way in which the probabilities Psoi(^'e ^e) and Pe,tm{i^n — ^i^) can fit 
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Summed Atmospheric Probabilities for Random Mass Matrices 




Figure 3: The sum of the probabihties Patm{j^n J^i) for atmospheric neu- 
trinos summed over i for i = e, fi,T for 10000 random mass matrices A4 are 
pfotted as a function of the parameter sinx^. 



these experimental results is if intcr-gcncrational mixing is suppressed so that 
z/g oscillates into rig and so that oscillates into n^. In other words, random 
mass matrices A4, even with sinxu ~ 0, produce probabilities Pgoii^e ^e) 
and Paxrai^ij. — ^^J) (suitably averaged respectively over the Earth's orbit and 
over the atmosphere) that are too small. The reason is that for small sinxy 
the flavor cigcnfield oi necessarily is split into two mass eigenfields, 
which reduces the probabilities to ~ |; and so, if there is further inter- 
generational mixing, then these probabihties tend to be too small. This 
effect is illustrated in Fig. 4 in which the probabilities Psoi(^e ~^ ^e) and 
Pa.tnl{^^l T^fj.) are plotted for 10000 ^ random mass matrices Ai all with 
sinXiy = 0.003. Figure 4 makes it clear that the probabilities Psoi(t'e ^e) 
and Patm(t'AJ J^n) are too small to fit the experimental results (66) and (67). 

The probabilities Psoi(^'e ^e) and Patm(^'/i ~^ tend to be somewhat 
larger when inter-generational mixing is limited. In Fig. 5 these probabilities 
are displayed for 10000 random mass matrices with smx,y = 0.003 and with 

^Due to the arXiv hmit of 650 kb on the length of articles, in some of the figures only 
6000 points are plotted. 
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Random Neutrino Mass Matrices with sin x^, = 0.003 



T 0.5 - 




0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 

^soKi^e Ve) 

Figure 4: The probabilities Psoi(i^e ^e) and -Patm(i^;u — i^n) for 10000 
random mass matrices Ai all with the parameter sin — 0.003. 



the singly off-diagonal matrix elements of D, E, and F suppressed by 0.2 and 
the doubly off-diagonal matrix elements suppressed by 0.004. The points in 
Fig. 5 are in much better agreement with (85) than are those of Fig. 4, but 
they still fail to match the experiments. Yet if the suppression of mixing is 
more severe, with factors respectively of 0.05 and 0.0025, then as shown in 
Fig. 6 the probabilities Pgoi(t'e — ^ ^e) and Patmi^fj. ^fi) do tend to cluster 
around (|, |) as required by the data. 

It is possible to relax the factors that suppress inter-generational mix- 
ing back to 0.2 and 0.04 and improve the agreement with the experimen- 
tal constraints (66) and (67) (while satisfying the CHOOZ constraint) pro- 
vided that one also requires that there be a quark-like mass hierarchy. The 
points in Fig. 7 were generated from 10000 random mass matrices M. with 
sinxi, — 0.003 and CKM-suppression factors of 0.2 and 0.04 as in Fig. 5 and 
with the i,j-th elements of the mass matrices E.F^ and D scaled by the 
factors f{i) * f{j) where / = (0.2, 1,2). Thus the mass matrix M. has the 
T, T elements that are larger than its /x, fi elements and /x, ii elements that 
in turn are larger than its e, e elements. The clustering of the probabilities 
-Psoi(i^e ^e) and Patm(i^/i — ^ i^//) around (|, |) in Fig. 7 shows that the ex- 
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Random Mass Matrices with sinaiy = 0.003 and Little CKM Mixing 
T 




0.3 0.35 0.4 



0.6 



Figure 5: The probabilities PgoK^^e ^e) and Pa.tra{^^i — for 10000 
random mass matrices Ai all with the parameter sinx,^ = 0.003 and with 
inter-generational mixing suppressed by factors of 0.2. 



Random Mass Matrices with sin.7;,y = 0.003 and Tiny CKM Mixing 
T 




0.35 0.4 0.45 



0.65 



Figure 6: The probabilities Psoi(^'e ^e) and Patmi^^i i^n) for 10000 
random mass matrices Ai all with the parameter sinXi, = 0.003 and with 
inter-generational mixing suppressed by factors of 0.05. 
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Random Quark-Like Mass Matrices with sin a;^ = 
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Figure 7: The probabilities Pso\{i^e J^e) and Patm(^'/i — i^n) for 10000 
random mass matrices Ai all with the parameter sinXj, = 0.003, with inter- 
generational mixing suppressed by factors of 0.2, and with a quark-like mass 
hierarchy. 

perimental results (66) and (67) are satisfied. The vector / was tuned so as 
to nearly saturate the cosmological upper bound (57) of about 8eV. 



Because is small, and because neutrinos oscillate mainly into sterile neu- 
trinos, the probabilities of the appearance of neutrinos are small, as shown 
by LSND and by KARMEN2. But since the mass differences among the 
three nearly degenerate pairs of neutrinos can lie in the eV range, such os- 
cillations should be observable. In Fig. 8 for a set of 10000 mass matrices 
generated randomly with the same parameters as for Fig. 7, the probabilities 
-P(i^M ^e) are shown for LSND and KARMEN2. 

In Fig. 9 for a set of 10000 mass matrices generated randomly with the 
same parameters as for Figs. 7 and 8, the probabilities P{i^f^ — > t'e) and 
-P(^/i i^fj.) for MiniBooNE are plotted. If intergenerational mixing is sup- 



LSND, KARMEN2, and MiniBooNE 
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LSND & KARMEN2 Probabilities for Random Matrices with sinx^ = 0.003 
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Figure 8: The probabilities PLSND(t'M ^ ^e) and PKARMEN2(i^/. ^ ^^e) for 
10000 random mass matrices Ai all with the parameter sinX;^ = 0.003, with 
inter- generational mixing suppressed by factors of 0.2, and with a quark-like 
mass hierarchy. 
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MiniBooNE Probabilities for Random Matrices with sinajy = 0.003 
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Figure 9: The probabilities PMiniBooNE(^'/. i^e) and PMimBooNE(t'M ^ J^fi) for 
10000 random mass matrices Ai all with the parameter sinx,^ = 0.003, with 
inter-generational mixing suppressed by factors of 0.2, and with a quark-like 
mass hierarchy. 



pressed only by factors of 0.2 and if MiniBooNE can achieve a sensitivity of 
0.001 for — > Pg and a precision of 0.01 for P^ — > P^, then it has a good 
chance of seeing both the appearance of Pg and the disappearance of P^. 

It may be, however, that intergenerational mixing is suppressed only by 
factors of 0.1. For as shown in Fig. 10, a set of 10000 matrices randomly 
generated, with CKM suppression factors of 0.1 for the off-diagonal terms 
and 0.01 for the off-off-diagonal ones and with the mass- hierarchy vector 
/ = (.2, 1, 2), agrees with the solar and atmospheric experiments just as well 
as those of Fig. 7. With these parameters the probabilities P(P// Pe) at 
LSND and KARMEN2 are smaller, as shown in Fig. 11. The chances of 
detection at MiniBooNE are similarly reduced as shown in Fig. 12. 
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Random Quark-Like Mass Matrices with sin = 0.003 
T 




0.65 



Figure 10: The probabilities Psoi(t'e — ^ ^'e) and Patm(^'fj ^'m) for 10000 
random mass matrices Ai all with the parameter sinx^ = 0.003, with inter- 
generational mixing suppressed by factors of 0.1, and with a quark- like mass 
hierarchy. 
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LSND & KARMEN2 Probabilities for Random Matrices with sinx^ = 0.003 
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Figure 11: The probabilities PLSNoi^f, ^ ^e) and PKARMEN2(i^,. ^ ^^e) for 
10000 random mass matrices Ai all with the parameter smx,y = 0.003, with 
inter-generational mixing suppressed by factors of 0.1, and with a quark-like 
mass hierarchy. 
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MiniBooNE Probabilities for Random Matrices with sina;^ = 0.003 
-T— q- 
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Figure 12: The probabilities PMimBooNE(^'M ^ ^e) and PMiniBooNE(t'/. ^ ^'^) for 
10000 random mass matrices Ai all with the parameter sinxj, = 0.003, with 
inter-generational mixing suppressed by factors of 0.1, and with a quark-like 
mass hierarchy. 
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Neutrinoless Double Beta Decay 



Neutrinoless double beta decay occurs when a right-handed antineutrino 
emitted in one decay n — p + e~ + z/g is absorbed as a left-handed neu- 
trino in another decay Ug + n ^ p + e^ . To lowest order these decays proceed 
via the Majorana mass term —iF*^ula2i'l. Let us introduce a second angle 
iji, defined by 



. 2 Tr(FtF) 

= Tr(£;t^ + FtF) - 

We have seen that we may fit the experimental data (66) and (67) by as- 
suming that smxi, ~ 0.003 and by requiring the mass matrices F, and D 
to exhibit quark-like mass hierarchies with little inter-generational mixing. 
Under these conditions the rate of Oi//3/3 decay is limited by the factor 

I Feel ^ ^ sin^ Xy siv? yy ml^ , (87) 

in which rrii,^ is the heavier of the lightest two neutrino masses. Thus the 
rate of Oi'f3f3 decay is suppressed by an extra factor ~ sin^ Xi, sin^ y^^ <, 10~^ 
resulting in lifetimes Ti q^^^ > 2x 10^'' yr. The B—L model therefore explains 
why neutrinoless double beta decay has not been seen and predicts that the 
current and upcoming experiments Heidelberg/Moscow, IGEX, GENIUS, 
and CUORE will not see Oi//3/3 decay. 



8 Conclusions 



The standard model slightly extended to include right-handed neutrino fields 
exactly conserves B ~ L if all Majorana mass terms vanish. It is therefore 
natural [8] to assume that the Majorana mass terms are small compared 
to the Dirac mass terms. An angle x^ is introduced that characterizes the 
relative importance of these two kinds of mass terms. When this parameter is 
very small, then the neutrinos are nearly Dirac and only slightly Majorana. 
In this case the six neutrino masses rrij coalesce into three pairs of nearly 
degenerate masses. Thus the very tiny mass differences seen in the solar 
and atmospheric neutrino experiments are simply explained by the natural 
assumption that x^, ^ or equivalently that B — L is almost conserved. 



33 



In these experiments the probabiUties ^801(^^6 ^ ^e) and Patm(^^t ^n) 
are respectively approximately one half and two thirds. One may fit these 
probabilities with random mass matrices in the eV range by setting sin x„ = 
0.003 and requiring the neutrino mass matrices E, F, and D to exhibit quark- 
like mass hierarchies with little inter-generational mixing. 

The three mass differences among the three nearly degenerate pairs are 
constrained only by the cosmological bound J2j ^ 8 eV and may lie in the 
range [0.01, 4] eV. 

This B — L model leads to these predictions: 

1. Because sin^a;,^ ~ and because inter-generational mixing is sup- 
pressed, neutrinos oscillate mainly into sterile neutrinos of the same 
flavor and not into neutrinos of other flavors. Hence rates for the ap- 
pearance of neutrinos, P{vi z/j/) with i 7^ i', are very low, as shown 
by LSND and by KARMEN2. But because the mass differences among 
the three nearly degenerate pairs of masses can lie in the eV range, 
MiniBooNE may detect ^ v^- 

2. The assumption that siv? Xy is very small naturally explains the very 
small differences of squared masses seen in the solar and atmospheric 
experiments without requiring that the neutrino masses themselves be 
very small. Thus the neutrinos may very well saturate the cosmological 
bound, fTT-j ~ 8 eV. In fact the masses associated with the points of 
Figs. 7-12 do nearly saturate this bound. Neutrinos thus may well be 
an important part of hot dark matter. 

3. The disappearance of u-,- should in principle be observable. 

4. In the B — L model, the rate of neutrinoless double beta decay is sup- 
pressed by an extra factor ~ sin^ Xi, sin^ y^, <, 10~^ resulting in lifetimes 
greater than 2 x 10^^ yr. Thus the current and upcoming experiments 
Heidelberg/Moscow, IGEX, GENIUS, and CUORE wiU not see Qiypp 
decay. 
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Appendix: LAPACK 

The thousands of 6 x 6 matrix computations displayed in the several graphs 
were performed in the inexpensive Lintel computing environment of Intel 
chips. Red Hat Linux, the Portland Group's Fortran compiler, and the linear- 
algebra software of the LAPACK collaboration. This appendix describes how 
to use LAPACK to perform the singular- value decomposition of an arbitrary 
6x6 complex matrix and the Takagi factorization of a 6x6 complex symmetric 
matrix. 

The driver subroutine ZGESVD [10] does a singular-value decomposition 
of an arbitrary matrix A. The FORTRAN 90 call is 

caU ZGESVD( JOBU, JOBVT, M, N, A, LDA, S, U, LDU, & 
VT, LDVT, WORK, LWORK, RWORK, INFO ) 

in which JOBU = 'A', JOBVT = 'A', M = N = LDA = LDVT = 6, WORK 
is an LWORK-dimensional double-complex vector, and RWORK is a 26- 

dimensional double-precision real vector. ZGESVD performs a singular- value 
decomposition of the 6x6 double-complex matrix A which is the mass matrix 
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M. and reports its singular values, which are the neutrino masses rrij, as 
the 6-dimensional double-precision real vector S. The matrices U and 
are contained in the 6x6 double-complex matrices U and VT. The integer 
INFO describes the level of success of the computation. The best value of the 
system-dependent integer LWORK is reported as the real part of WORK (1); 
for the computations of this paper, which were done on Pentium II's and Ill's, 
LWORK was 396. If JOBU is 'A' and JOBVT is 'A', then the subroutine 
call destroys the matrix A. Thus one must redefine A before every call to 
ZGESVD. 

If the matrix A is symmetric, then one may convert the singular-value 
decomposition A = U SV"^ into the Takagi factorization A = WSW^ by the 
FORTRAN 90 equivalent of 



W, 



i 



^jk \Ukj\ 



(88) 



where the index k is chosen to avoid any possible singularity. 
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